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Abstract. In this article we study the sub-Riemannian geometry of the 
spheres S 2n+1 and S 4n+3 , arising from the principal S 1 — bundle structure 
defined by the Hopf map and the principal S 3 — bundle structure given 
by the quaternionic Hopf map respectively. The S 1 action leads to the 
classical contact geometry of S 2n+1 , while the S 3 action gives another 
type of sub-Riemannian structure, with a distribution of corank 3. In 
both cases the metric is given as the restriction of the usual Riemannian 
metric on the respective horizontal distributions. For the contact S 7 
case, we give an explicit form of the intrinsic sub-Laplacian and obtain 
a commutation relation between the sub-Riemannian heat operator and 
the heat operator in the vertical direction. 

1. Introduction 

One of the main objectives of classical sub-Riemannian geometry is to 
study manifolds which are path-connected by curves admissible in a certain 
sense. Admissibility refers to a constraint on the velocity vector of an ab- 
solutely continuous curve 7 : [0, 1] — > M, where M is a smooth connected 
manifold. More precisely, if H G TM is a smooth distribution, then 7 is 
admissible or horizontal if j(t) G % a.e. The distribution H is often called 
horizontal distribution in the literature. 

The idea of studying sub-Riemannian geometry arising from well-behaved 
fiber bundles was introduced by R. Montgomery in [T5], although the Rie- 
mannian analogue had been studied many decades before. The idea is the 
following: given a submersion tt : Q — > M between two Riemannian man- 
ifolds Q and M, where dimM < dimQ, define a "horizontal" distribution 
over Q by the pull-back bundle ir*(TM) of the tangent bundle of M via tt. 
In the case when we have a principal G— action over Q preserving the fibers 
of the submersion, the manifold M can be identified with the orbits of the 
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action and, after some technical assumptions, it is possible to obtain an 
explicit characterization of sub-Riemannian geodesies. 

The aim of the present article is to describe the sub-Riemannian geom- 
etry of two sub-Riemannian structures for odd-dimensional spheres. More 
specifically, we study the sub-Riemannian geometry arising from the contact 
distribution for the spheres S 2n+1 with metric given as a restriction of the 
usual Riemannian metric, and the one arising from the quaternionic Hopf 
fibration for the spheres S 4n+3 . 

This article is organized as follows. In Section 2, we give some standard 
definitions of sub-Riemannian geometry which will be needed in the rest of 
the paper. In Section 3 we give an explicit description of sub-Riemannian 
geodesies in spheres S 2n+1 endowed with the standard contact distribution 
and we study some of their geometric properties. In Section 4 we use the 
obtained form of geodesies in the case of S 3 to give another interpretation to 
a result by Hurtado and Rosales in [TT] . With this new point of view, we are 
able to extend their result to contact spheres of an arbitrary odd dimension. 
Section 5 is the analogue to Sections 3 and 4 for the case of spheres of the 
form S 4n+3 endowed with a distribution of corank 3. Section 6 is somewhat 
different technically, but it is in spirit related to the core of this article. 
It deals with a geodesic differential equation for the quaternionic H— type 
group studied in [4] , obtained generalizing the techniques in [17] . The reason 
for studying this equation here is to pose the question of a similar equation 
for the case of S 7 and a distribution of rank 4. Section 7 consists of the 
construction of the intrinsic sub-Laplacian for S 7 . The main result states 
that it is the sum of the squares of an orthonormal basis of the horizontal 
distribution. Finally, Section 8 employs the previous construction to obtain 
a simple form of the heat operator for S 7 in a similar way as obtained in [3] . 



2. Preliminaries and notations 

2.1. Sub-Riemannian geometry. Let us first give some general defini- 
tions, which will be adapted to our purposes when it will be necessary. Let 
M be a smooth connected manifold of dimension n, together with a smooth 
distribution H C TM of rank k, 2 < k < n. The manifolds of our interest 
are endowed with distributions satisfying the bracket generating condition, 
i.e. distributions whose Lie hull equals the full tangent bundle of M. To be 
more precise, define inductively the vector bundles 

H 1 = H, H r+1 = [H r , H] + W for r > 1, 

which naturally induce the flag 
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We say that H is bracket generating if for all x G M there is an r(x) G Z 
such that 



If the dimensions dim?^. do not depend on x for any r > 1, we say that H 
is a regular distribution. The least r such that (JTJ is satisfied is called the 
step of %. In this paper we will focus on regular distributions of step 2. 

A natural question to pose is, given M and H, whether one can join any 
two points of M via a horizontal curve, i.e. an absolutely continuous curve 
7 : [0, 1] — >■ M which satisfies j(t) G H almost everywhere. A complete 
answer to this question was given in [19], which shows a deep generalization 
the celebrated Chow-Rashevskii theorem, see [3116], that gives a sufficient 
condition and can be stated as follows: 

Theorem 1. Let M be a connected manifold and H C TM be a bracket 
generating distribution, then the set of points that can be connected top G M 
by a horizontal path coincides with M . 

Remark: A slightly more general version of Theorem [1] states that, if M is 
not connected, then the set of points that can be connected to p G M by a 
horizontal path is the connected component containing p. Since we assumed 
the manifold to be connected, the general formulation is unnecessary. 

After these preliminaries, we are ready to specify the class of manifolds 
of our interest. 

Definition 1. A sub-Riemannian structure over a manifold M is a pair 
(1-L,(-,-)sr)> where H is a bracket generating distribution and (-,-) s r is a 
fiber inner product defined on 7-L. The triple (M,7i, (■,-) s r) is called sub- 
Riemannian manifold. 

In this context, the length of a horizontal curve 7 : [0, 1] — > M is defined 
to be 



where || 7(£)|| 2 = (j{t), j(t)) S R whenever ■jit) exists. 

This notion of length gives rise to the Carnot-Caratheodory distance 
d(p,q) between two points p,q G M, given by d(p,q) := mi£(j), where 
the infimum is taken over all absolutely continuous horizontal curves joining 
p to q. An absolutely continuous horizontal curve that realizes the distance 
between two points is called a horizontal length minimizer. It is clear that 
if % is bracket generating then d(p, q) is a finite nonnegative number. 

Considering a trivializing neighborhood U p around p G M for the sub- 
bundle "H, one can find a local orthonormal basis Xx, . . . ,Xf. with respect 




H r x (x) = T X M. 
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to (•, -) s r- The associated sub-Riemannian Hamiltonian is given by 

k 

H{qA) = \Y J KX m (q))\ 

m=l 

where (q, A) G T*U P . A normal geodesic corresponds to the projection to 
Up C M of the solution of the Hamiltonian system 

dH 

o\i 

■ OH 

Xi — q — , 

dqi 

where (q i7 Xi) are the coordinates in the cotangent bundle of M. 
Remark: It is possible to define sub-Riemannian geodesies in a more general 
context. There are many interesting problems related to the classification 
of such curves, their analytic and geometric properties. In [12] the prob- 
lem for the case of rank two distributions is studied and essentially solved. 
Nevertheless, in the case of step two distributions, the general notion of 
geodesic gives rise to two cases: curves consisting of one point and normal 
geodesies. Thus, normal geodesies are the only interesting case for our pur- 
poses. Note that in this case normal geodesies are local length minimizers, 
in the sense that any sufficiently small arc of a normal geodesic minimizes 
the length functional. On the other hand one of the particular features of 
sub-Riemannian geometry, as the sub-Riemannian Heisenberg group exem- 
plifies, is that it is possible to find arbitrarily close points that can be joined 
by normal geodesies with different lengths. 

2.2. Sub-Riemannian principal bundles. Our first goal is to recall a full 
characterization of normal geodesies in the case of sub-Riemannian princi- 
pal bundles. As a direct application we obtain an explicit formula for the 
sub-Riemannian geodesies on odd-dimensional spheres, with respect to dis- 
tributions of corank 1 and 3 in Sections 3 and 5 respectively. For the sake 
of completeness we recall some definitions and notations given in [T5] . 

For a submersion tt : Q — > M with fiber Q m = tt (m) through m € M, 
the vertical space at q G Q is given by T q Q w ^ and it is denoted by V g . In 
this context, an Ehresmann connection for 7r : Q — > M is a distribution 
H C TQ which is everywhere transversal to the vertical space, that is: 

V q ®H q = T q Q for every q E Q. 

Let us assume that a Lie group G acts on Q in such a way that it : Q — > M 
becomes a fiber bundle with fiber G. We say that the submersion it is a 
principal G— bundle with connection % if the following conditions hold: G 
acts freely and transitively on each fiber, the group orbits are the fibers of 
7i : Q — > M, and the G— action on Q preserves the connection "H. Observe 
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that the second condition implies that M is isomorphic to Q/G and 7r is the 
canonical projection. We will refer to the connection % as the horizontal 
distribution. 

For the rest of this section, let us denote the Lie algebra of G by g, and 
the corresponding exponential map by exp G : g — > G. 

Definition 2. For the principal G— bundle n : Q — )■ M, the infinitesimal 
generator for the group action is the map a q : g — > T q Q defined by 



».<f) - i 



?exp G (eO 



for q G Q and £ G g. If the metric (■, •) in Q is G— invariant, we have a 
well-defined bilinear form 

which is called the moment of inertia tensor at q. 

The G— invariant Riemannian metric on Q is said to be of constant bi- 
invariant type if its moment of inertia tensor I q is independent of q G Q. 
Recall also that, in the case of a principal G— bundle, for each q G Q the 
infinitesimal generator a q is an isomorphism between the vertical space V q 
and g. We refer to its inverse as the g valued connection one form. 

With all of these at hand, we can state the main tool required in this 
section. This will imply almost immediately Corollaries Q] and [2] which are 
of core importance in the present paper. The proof of the following theorem 
can be found in IT51. 



Theorem 2 (Horizontal Geodesies for Principal Bundles) . Let n : Q — )■ M 

be a principal G— bundle with a Riemannian metric of constant bi-invariant 
type. Let 7-L be the induced connection, with g valued connection one form A. 
Let exp R be the Riemannian exponential map, so that 7a(t) = exp R (tv) is 
the Riemannian geodesic through q with velocity vector v G T q Q. Then any 
horizontal liftj of the projection tto'j r is a normal sub-Riemannian geodesic 
and is given by 

j(t) = exp R (tv) exp G (-tA(v)) 
where exp G : g — > G is the exponential map of G. Moreover, all normal 
sub-Riemannian geodesies can be obtained in this way. 

Remark: In Theorem [21 the sub-Riemannian geodesies are considered with 
respect to the metric induced by restricting (•, •) to %. Recall that constant 
bi-invariant metrics must be G— invariant. 

3. Sub- Riemannian Geodesics on S 2n+l 

In the case of odd dimensional spheres S 2n+1 , embedded as the boundary 
of the unit ball in C™ +1 , there is a natural action of S 1 = SU(1) on it, 
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via componentwise multiplication by a complex number of norm 1. This 
action induces the well known Hopf fibration S 1 — > S 2n+1 — > CP™, which 
forms a principal S 1 — bundle with connection H given by the orthogonal 
complement to the vector field 

(2) V n+1 (p) = -y d X0 + x d yo - ... - y n d Xn + x n d Vn 

at each p = (x , yo, ■ ■ ■ , x n , y n ) G S 2n+1 , with respect to the usual Riemann- 
ian metric of S 2n+1 as embedded in M 2 (™ +1 ) S C n+1 . In [9] it is shown that 
this distribution coincides with the holomorphic tangent space HS 2n+1 of 
S 2n+l thought as an embedded CR manifold and that it also coincides with 
the contact distribution given by keru; with respect to the contact form 

u = -y dx + x dy - ... - y n dx n + x n dy n . 

Note that the components of the vector V n+ i(p) are the same as in the su(l) 
action i ■ p. 

As a direct application of Theorem El it is possible to describe all sub- 
Riemannian geodesies for the sphere S 2n+1 as a sub-Riemannian manifold 
equipped with connection 7i and with metric restricted from M 2 ( n+1 ). By 
the results discussed in [5], the holomorphic tangent space for S 2n+1 is the 
distribution induced by the principal S 1 — bundle given by the Hopf fibration 
S 1 — > S 2n+l — > CP™ with su(l)— valued connection form A{v) = i(v, V n+ i), 
v G T p S 2n+1 , V n+ \ denotes V n+ i(p) and (•, •) stands for the standard inner 
product in IR 2(n+1 ). Moreover, the usual Riemannian structure on S 2n+1 is 
of constant bi-invariant type, since we have 



d 
de 



?exp su(1) (eO = ai ■ q = aV n+1 (q), 



(i) 

e=0 



for any q G S 2n+1 and £ = ia G su(l). Therefore, the inertia tensor is given 

by 

I q (ia,ia) = (aV n+1 (q) , aV n+1 (q)) = aa, 

which does not depend of the point. 

By Theorem El we have the following result. 

Corollary 1. Let p G S 2n+1 = {(z , . . . , z n ) G C n+1 : \z \ 2 + . . . + \z n \ 2 = 1} 
and v G T p S 2n+l . If 7h(£) = (zo(t) , . . . , z n (t)) is the great circle satisfying 
7r(0) — P an d 7_r(0) = v , then the corresponding sub-Riemannian geodesic 
is given by 

(3) j(t) = {z (t)e- u ^ v ^\ . . . , z n (t)e- u ^ v ^) . 

In order to analyze in more details formula (j3J), let us introduce some 
notations and the necessary setup. Recall that the Riemannian geodesic 
starting at p G S n with velocity v G T p S n of any sphere S n as a submanifold 



SUB- RIEMANNIAN ODD DIMENSIONAL SPHERES 



7 



of M n+1 , with the standard Riemannian structure, is given by: 

(4) lR (t) =pcos(\\v\\t) + -^- S m(\\v\\t), 

where \\v\\ 2 = (v,v). In the case of our interest, a great circle j R (t) in S 2n+l 
as a submanifold of IR 2(n+1 ) = C n+1 will be written in complex notation as 
7r(£) = ( z o(t) , . . . , z n (t)) . For notational simplicity, the action of A G S 1 
over (p , . . . ,p n ) G S 2n+1 is denoted by A • p = (Xpo, • • • , Ap n ). Let us write 
7(0) = 7i?(0) = p = (ao + i&o, ..-A + ^n) G S 2n+1 and 7a (0) = u = 
Oo + 0o, ...,«„ + i/3 n ) e TpS 12 " 4 " 1 . Observe that K,+i(7(£)) = i • t(*)- As 
above, V^+i = V^+i( 7 (0)). 

Remark: In the subsequent calculations, the notation (•,•)# will denote 
the standard Hermitian product in C™ +1 . We recall that the standard inner 
product (•, •) in R 2 (™ +1 ) satisfies 

Re (;-) H =(;■)■ 

Theorem [2] assures that 7 is a horizontal curve, i.e. (j(t), V^+i( 7 (t))) = 0, 
nevertheless it is possible to check directly this by straightforward calcula- 
tions. Since some of the computations will appear later, it is convenient to 
write them down. First notice that 

(7(*),W7(*)))h = ((^(v,V n+1 ) lR (t)+^ R (t))e-^ v ^\ 

ie -^y^ lR it)) H 

= -{v, V n+1 ) ( f y R (t),j R (t)) H -i(jR(t),j R (t)) H 
= -(v, V n+1 ) - i{j R (t),j R (t)) H . 

Thus the problem is now to determine the value of 

n 

( A f R {t), lR {t)) H = j2 i k{t)'^t). 

k=0 

By straightforward calculations, it is easy to see that 

n n 

^z k (t)zi$j = (cos 2 (\\v\\t)-sm 2 (\\v\\t))J2( a kak + hf3 k ) + 

k=0 k=0 

n 

+i^2(a k (3 k - b k a k ) 

k=0 

= {p,v)cos(2\\v\\t)+i(v,V n+1 ) 

(5) = i(v,V n+1 ), 

yielding to (j(t), V^ + i(7(t)))ij = 0, which implies the horizontality of the 
curve j(t). 

Let us now address the problem of connecting two points in S 2n+1 by sub- 
Riemannian geodesies. We know by Theorem [T] that it is possible to find a 
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horizontal curve T : [0,T] — > S 2n+l such that 

(6) T(0) = p and T(T) = q, 

for any pair p,q G S 2n+1 and all fixed time parameter T > 0. A natural 
question to ask is whether T can be taken as a geodesic in (jf)]). Due to the 
complexity of the problem, we will give a partial answer to it. It is important 
to remark that Proposition [T] is a direct analogue of the result obtained in 
(5j Theorem 1] in the particular case of n = 1, i.e. for the three dimensional 
sphere. 

Proposition 1. The set of sub-Riemannian geodesies arising from great 
circles ■jnit) such that 7r(0) G % = kerw is diffeomorphic to CP n . 

Proof. In this case any sub-Riemannian geodesic starting at p G 5' 2n + 1 with 
initial velocity v G % C T p S 2n+1 coincides with the corresponding great 
circle, since the condition 7r(0) G % — kero; is equivalent to (v, V n+ i) = 0, 
thus 

T (f) =pcos(Ht) + ^sm(|M|t) 

whose loci is uniquely determined by the point [v] G CP n . □ 

Observe that this CP™ can be seen as a submanifold of S* 2n+1 which is 
transversal to V n+ \ along the fiber containing p. As remarked in [5] for S 3 , 
this can be seen as a sophisticated analogue of the horizontal space at the 
identity in the (2n + 1)— dimensional Heisenberg group. 

Let us conclude this discussion with an interesting result which will be 
of importance in the following Section. This can be thought of as a sort of 
Pythagoras theorem for contact spheres. 

Proposition 2. For a horizontal sub-Riemannian geodesic of the form 

7 (f) = {z (t)e- u ^ v ^\ . . . , z n (t)e- u ^ v ^) 
the following equation holds 

|| 7 (f)|| 2 + (v,V n+1 ) 2 =\\v\\ 2 . 
Thus, its velocity is constant and its sub-Riemannian length for t G [a, b] is 
£(7) = (6-a)VIH| 2 -(^K + i) 2 - 
Proof. By straightforward calculations, we have 

<7(*),7(*)>h = ((-i(v,V n+1 ) lR (t)+Ut))e- l{vy " +l)t , 
(-i(v, V n+1 ) lR (t) + 7«(f))e-''^' y " +l)i )H 
= {v,V n+1 ) 2 ( lR (t), lR (t)) H + (7*(*),7ji(*)>h 

+{v,V n+ i)(i('j R ,'y R )H - i{lR,iR)H) 
= (v,V n+1 ) 2 + \\v\\ 2 -2(v,V n+1 ) 2 . 

Here we have used equation (JSJ). The proposition follows. □ 
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Remark: According to Proposition |2j the condition that a curve j(t) = 
e~ tt ^ v ' Vn+1 ' r yji(t) is parameterized by arclength is equivalent to require that 
IMI 2 = l + (*;, K+i> 2 . 

4. Curvature of sub-Riemannian geodesics on S 3 

In [TT], the authors describe the horizontal geodesics of the three dimen- 
sional sphere with respect to its contact distribution, obtaining an explicit 
expression for these curves. The key tool to achieve this is the following 
proposition. 

Proposition 3. Let 7 : / — > S 3 be a C 2 horizontal curve parameterized by 
arc-length. Then 7 is a critical point of length for any admissible variation 
if and only if there is A G M such that 7 satisfies the second order ordinary 
differential equation 

(7) V^ 7 + 2A J( 7 ) = 0, 

where V is the Levi-Civita connection and J is the standard almost complex 
structure on S 3 . 

The authors call the parameter A above the curvature of 7, since after 
projecting it via the Hopf fibration, A becomes precisely the curvature of the 
projected curve in S 2 . Note that the curves with zero curvature are precisely 
the horizontal great circles. It is our purpose to find an explicit expression 
for A in terms of known parameters of the sub-Riemannian geodesics of S 3 , 
as presented in Corollary [TJ 

Proposition 4. The curvature of the sub-Riemannian geodesic 

l{t) = e-^ l lR {t) 
in S 3 , parameterized by arc-length, equals (v, V 2 ). 

Proof. The Lie group structure of S 3 as the set of unit quaternions, induces 
the globally defined vector fields 

V(p) = -Vid Xl + x t d yi - y 2 d X2 + x 2 d y2 , 

(8) X(p) = -x 2 d Xl + y 2 d yi + x x d X2 - yid V2 , 
Y{p) = ~V2d Xl - x 2 d yi + y x d X2 + x x d V2 , 

at p = (xi,yi, x 2 , y 2 ) G S 3 , which are orthonormal with respect to the usual 
Riemannian structure of 1R 3 . Observe that V(p) = V 2 (p) as defined in (j2J). 

Let p = (xi,yi,x 2 ,y 2 ) = 7(0) G S* 3 be the initial point of 7 and let v = 
(v Xl ,v yi ,v X2 , v y2 ) = 7r(0) G T P S 3 be the initial velocity of the corresponding 
great circle. By direct calculation, we have 

(9) j(t) = fx(t)X(j(t)) + f Y (t)Y(j(t)), 
where, denoting a = (v, X), (3 = (v, Y), we have 

f x {t) = a cos(2t(t», V)) + p sin(2t(t», V)), 
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f Y {t) = (3 cos(2t(v, V)) - a sin(2t(v, V)). 
It follows from this decomposition that 

(10) J(7(t)) = -f Y (t)XMt)) + fx(t)Y(j(t)). 

It remains to determine the term Vyy. It is well-known that for sub- 
manifolds of R™, the vector field V^,7 corresponds to the projection of the 
second derivative 7 to the tangent space of the submanifold. In this case, 
differentiating we obtain 

V^ 7 = 2(v,V)(f Y (t)X( 1 (t))-f x (t)Y( 1 (t))) 

= -2(v,V) J{>y(t)). 

The proposition follows. □ 

Remark: Note that in case p = (1, 0, 0, 0) G S 3 , a great circle starting at p 
with velocity vector v = (0,v yi ,v X2 ,v V2 ) G T P S 3 is given by 

lR {t) = (cos(\\v\\t),^-sm(\\v\\t), ^sm(|M|t), ^sm(\\v\ 

V IMI \\ v \\ \\ v \\ 

Then, the corresponding sub-Riemannian geodesic is 

(11) 7 (t) = e-^nit), 

where v^+v^ = 1, since the curve is parameterized by arc- length. It follows 
that the curvature is given by (v, V 2 ) = v yi . 

In [TT] the problem of existence of closed sub-Riemannian geodesies is also 
discussed. Their result is that a complete geodesic 7 in S 3 parameterized 
by arc- length, with curvature A is closed if and only if A/a/1 + A 2 G Q. This 
result can be generalized to any odd dimensional sphere. 

Proposition 5. Let 7 : IR — > S 2n+1 be a complete sub-Riemannian geodesic 
parameterized by arc-length, with initial velocity v G T p S 2n+1 . Then 7 is 
closed if and only if 

- = G Q 

y/l + (v,V n+1 )* 

Proof. The curve 7:^—7' S 2n+1 is closed if and only if for some T > 
i(v,v n+1 )T f pcOB q v ^ + JL s i n ^\ v \\ T ) 



p = e 

bmce v G T p S 2n+ \ we know that v is orthogonal to the vector joining 
G M. 2n+2 to p, with respect to the usual Riemannian structure of R 2n+2 . 
This means that sin(||t> ||T) = 0, which forces T = kir/\\v ||, k G Z. 

To complete the argument, we only need to see that 

± e -ik({v,V„+i)/\\v\\yn-p _ p 
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if and only if 

(V, V n+l) = (v,Vn+l) e q 

\\v\\ y/i + (v,v n+1 y ' 

where we have used the remark after Proposition [2J □ 

5. SUB-RlEMANNIAN GEODESICS ON S 4n+3 

Let us consider the sphere 5" 4n+3 embedded as the boundary of the unit 
ball in (n + 1)— dimensional quaternionic space IEP +1 . As usual, let us denote 
the quaternionic units as i, j, and k. There is a natural right action of 
Sp(l) = S 3 on H n+1 , via componentwise multiplication by a quaternion of 
norm one. This action induces a quaternionic Hopf fibrations S 3 — > S 4n+3 — > 
HLP n , given by 

H : S 4n+3 ->■ HP n 

(g ,-..,g„) ^ [q : . . . : q n \. 

This submersion forms a principal S 3 — bundle with connection given by 
the orthogonal complement to the vector fields 

v n+i{p) = -yod Xo +x dy +w d ZQ -z d Wo -. . .-y n d Xn +x n dy n +w n d Zn -z n d Wn , 
V u+i(p) = -z d xo -w d yo +x d zo +yod wo -. . .-z n d Xn -w n d yn +x n d Zn +y n d WnJ 
v n+i(p) = -wod Xo +z dy -yod Zo +x d Wo -. . .-w n d Xn -z n d yn +y n d Zn +x n d Wn , 

at each p = (xo,yo, Zq,Wq . . . ,x n ,y n , z n ,w n ) G 5' 4n+3 , with respect to the 
usual Riemannian metric of S in+3 as embedded in R 4(n+1 ) = H ra+1 . It is easy 
to see that the following commutation relations hold for V^ +1 , V^ +1 , V 3 +l 

[t&i, K+i] = 2K+D [K+i, v 3 +1 ] = 2V n x +1 , [V* +1 , V 3 +1 ] = -2V* +1 . 

Thus one recovers the fact that spa.n{V^ +1 (p) , V^ +1 (p), V 3 +1 (p)} is isomor- 
phic as Lie algebra to sp(l), the Lie algebra associated to S 3 . 

It is a well established fact that this distribution is bracket generating. 
In fact, the geometry of this spheres S An+3 is known to be a quaternionic 
analogue of CR-geometry, see [2]. Note that the components of the vector 
V^ +1 {p) are the same as in the sp(l) action p ■ i. Similar statements hold for 
V^ +1 (p), V 3 +l {p) and p ■ j, p ■ k respectively. 

In order to apply Theorem [2] in this situation, it is necessary to specify 
the sp(l)— valued connection form associated to the submersion H. In this 
case, the connection form is given by 

A(v) = i{v, Ci) + j(v, K+i> + k{v, V 3 +1 ) 

where v G T p S 2n+1 , V" +1 denotes V^ +1 {p) (a = 1,2,3) and (•, •) stands for 
the standard inner product in M 4(n+1 ). Moreover, the usual Riemannian 
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structure on S 4n+3 is of constant bi-invariant type, since for any q G S An+3 
and £ = za + j(3 + G sp(l), a, /3, 7 G K we have 

-7- 9 ex P« P (i)(^) = aq-i + pq-j + wk 

Therefore, the inertia tensor is given by 

I q (ia + j/3 + &7, ia + j/3 + £7) = 

= (aK +1 (^^ +1 (g) + 7 ^ +1 ( ? ), 5K +1 ( 9 )^ +1 (g) + ^(g)) = 

= act + f3(3 + 77, 

which does not depend of the point. 

As for Corollary [TJ we have the following result. 

Corollary 2. Letp G S in+3 = {(u ,...,u n ) G H n+1 : \u \ 2 + . . . + \u n \ 2 = 1} 
and v G T p S An+3 . If jn(t) = (u (t), . . . ,u n (t)) is the great circle satisfying 
7i?(0) = p and 7r(0) = v , then the corresponding sub-Riemannian geodesic 
is given by 

(12) 7 (f) = (««,(*) • e-^W • • .,«„(*) • e- u <»)) . 

In Corollary [21 the quaternionic exponential is defined by 

e a i+bj+ck = v / fl2 + 6 2 + C 2 + gin v / fl 2 + 6 2 + C 2 . ^ + ^ + Cfc 



Va 2 + 6 2 + c 2 

for a, 6, c G K. Note that the curve e~ tA ^ is simply the Riemannian geodesic 
in S" 3 starting at the identity of the group e = (1, 0, 0, 0), with initial velocity 
vector (0, -(v, V^i), V% +1 ), -(v, V 3 +l )). 

Corollary [5] implies immediate analogues to Proposition [T] and to Proposi- 
tion 0, which we state for the sake of completeness. Proofs are adaptations 
of the aforementioned Propositions. 

Proposition 6. The set of sub-Riemannian geodesies in S 4n+3 arising from 
great circles 7_r(£) such that 7r(0) is orthogonal to V£ +1 , V 2 +1 and V 3 +1 is 
dijfeomorphic to HP". 

Proposition 7. Let 7 : IR — >■ S 4n+3 be a complete sub-Riemannian geodesic 
parameterized by arc-length, with initial velocity v G T p S 2n+1 . Then 7 is 
closed if and only if 

(y, YM YM (<YM 

II 119 ' II 119 ' II 119 
\\ V \\ \\ V \\ \\ V \\ 

In analogy with Proposition |2} let us consider a similar statement in the 
case of the spheres S 4n+3 . 
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Proposition 8. For a horizontal sub-Riemannian geodesic of the form 

1 (t) = {w (t)-e- tA ^,...,w n (t)-e- tA ^) 
the following equation holds 

\\i(t)\\ 2 +\\A(v)r = \\v\\ 2 , 

where \\A(v) || 2 = (v, V^f + (v, V' +1 f + (v, V 3 +l )\ 

Proof. Recall that if 7 is a sub-Riemannian geodesic, then the length of the 
velocity vector ||7(t)|| does not depend on t. Thus without loss of generality 
we can assume t = 0. Let us introduce the following notation 

p = -f(0) = (x ,y ,z ,w ,...,x n ,y n ,z n ,w n )eS in+3 , 

v = 7i?(0) = (v Xo ,v yo ,v ZQ} v Wo} . . . , v Xn ,v yn ,v Zn ,v Wn ) G T p S* n+3 . 

Differentiating equation ( |T2|) and evaluating at t = 0, we have 

7(0) = v-(v, KVi^Vi - («, - (v, v 3 +1 )v 3 +1 . 

The orthogonality of the vector fields V l +1 , V^ +1 , V^ +1 implies the desired 
relation. 

□ 

6. Curvature of sub-Riemannian geodesics on H 1 

The proof of Proposition [3] is given in [T7j for the case of the three di- 
mensional Heisenberg group. As mentioned in [TT], the proof for the case 
of the sub-Riemannian three dimensional sphere is basically the same. The 
authors have pointed out, in private communication, that the same result 
holds for all three dimensional pseudo-Hermitian manifolds. 

Note that that if M is either the Heisenberg group of topological dimension 
3 or the sphere S 3 , with Reeb vector field R, then the quotient vector bundle 

TM/span{R} M 

is trivial. We have not been able to show that the corresponding vector 
bundle 

TS 7 /s pa n{V 2 \V 2 2 ,V 3 }^S 7 

is trivial, which makes difficult to find an analogous argument to the one 
employed in [TTj . 

The main goal of this section is to find an analogue to Proposition [3] for the 
Gromov-Margulis-Mitchell-Mostow tangent cone of S 7 , see [101 [131 [TH [15] . 
which corresponds to the seven dimensional quaternionic if— type group H 1 , 
as presented in [1]. Observe that the idea of studying the tangent cone before 
the sub-Riemannian manifold of interest corresponds to the case in [TTj, 
since the three dimensional Heisenberg group is the tangent cone to the 
sub-Riemannian S 3 . We will study whether this method extends to S 7 in a 
forthcoming paper. 
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6.1. The quaternionic H— type group H 1 . Let us consider the 4x4 
matrices X, J and /C, given by 



/ 





1 





\ 






( 








-1 \ 




-1 











, J 










-1 
















1 







1 








\ 





- 


-1 


°) 






V 1 








0/ 










(0 


- 


-1 


\ 

























1 












K 




1 

























v° 


-1 














Note that X, J and /C are a fixed representation of the quaternion units, i.e. 
if U denotes the identity matrix of size 4x4, then span{W,X, J~, K.} = HI as 
algebras via the isomorphism 

up : span{W,X, J,JC} ->■ H 

given by <p(U) = 1, <£>(X) = i, y?(j7") = j, <p{fc) — k and extended by linearity. 

The seven dimensional quaternionic if— type group H 1 corresponds to the 
manifold M 4 © M 3 with the group operation o defined by 

(x, z) o {x\ z') = + x', zx + z' x + - x' T Xx, 

Z J + z J + \ x ' T J x i zk + z'k + t^ x' t K,x^ 

where x, y, z are column vectors and x' T , y' T , z' T are row vectors in IR 4 . 

The Lie algebra f) 1 corresponding to H 1 is spanned by the left invariant 
vector fields 



Xi(x,z) 
X 2 (x,z) 
X 3 (x,z) 
X 4 (x, z) 
Z x {x, z) = 



d 1 
dx^ + 2 

d 1 

dx 2 ' + 2 

d 1 
dx~ 3 + 2 

d 1 

dx 4 2 
9 



+x 2 



-Xi 



+X4 



-X 3 



d 

dz T 
d 

dz x 
d 

dz x 
d 
dz T 



X4 



- x 3 



+ x 2 



Xi 



d 

dzj 

d 

dzj 

d 
dzj 
d 



X3 



+ x A 



+ Xi 



d 

dzic 

d 
dz K 

d 
dz K 

d 



%2- 



Zj(x,z) 



d 



dzj dz K 
Z K (x,z) 



d 



ozx dzj dz K 

at a point (x,z) = (x±, x 2 , x 3 , x 4 , zx, zj, zjc) G H 1 . A Riemannian metric 
(•, •) in H 1 is declared so that X ll . . . , X 4 , Z X) ■ ■ ■ , Z^ is an orthonormal frame 
at each (x, z) e H 1 . The sub- Riemannian structure on H 1 we are interested 
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in is denned by the left invariant distribution T> = spanjXx, X 2 , X 3 , X 4 } and 
the restriction of the metric previously defined. 

Observe that T> is bracket generating of step two. In fact, we have the 
commutator relations 

[Xi,X 2 ] = [X 3 ,X 4 ] = — 

(13) [X 2 ,X 3 ] = [X 1? X 4 ] = Zj, 

[X 1 ,X 3 \ = [X 4 ,X 2 ] = Z K . 

All the remaining commutators between the chosen basis of f) 1 vanish. 

From the well-known Koszul formula for the Levi-Civita connection asso- 
ciated to the metric (-, •) 

(Z, VyX) = \{X{Y, Z) + Y(Z, X) - Z(X, Y)- 

-([X, Z], Y) — ([Y, Z], X) - ([X, Y],Z)), 

see for example [8], the orthonormality of the basis {Xi, . . . , X 4 , Zx, ■ ■ ■ , Z/c}, 
and equations f lT3|) we get that 

(X b ,V Xa Z r ) = -l([X a ,X b ],Z r ), (Z s ,V Xa Z r ) = 0, 

for any a,b = 1, . . . , 4, r, s = X, J , K. This translates to the equation 

1 4 

(14) V Xa Z r = --J2([Xa,X b ],Z r )X b , 

b=i 

which reduces to the following identities 
Vxi^x = -X 2 , Vx 2 Zi = _ 2^' ^ x z Z x = 2"^ 4 ' ^x 4 Zx = ~2^ 3 ' 

^ 1 x x Zj = -t^ 4 ' ^x 2 Zj = _ 2^ 3 ' ^ x 3 Zj = 2^2, ^ x A Zj = -Xi, 

Vxi^/c = — 2^ 3 ' ^ x 2 Zk, = 2^ 4 ' ^x 3 Zk. = 2^ 1 ' ^x 4 Zic = ~2^ 2 ' 
Therefore, it follows that the maps J r : T> — > T> defined by 
J r (X) = 2VxZ r , r = l,J,)C, 
are almost complex structures. Note that the equation 

(15) (Jr(U 1 ),U 2 ) + (U 1 ,J r (U 2 )) = 

holds for every r = X, J7", K, and every Ui,U 2 G T>. Note in particular that 
equation (TIoT) implies that (U, J r {U)) = for all U G V. 
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6.2. A variational argument. Consider a manifold M and let H C TM 
be a distribution. A variation of a curve 7 : [a, b] — > M is a C 2 -map 7 : 

x I 2 — > M, where 1%, J 2 are open intervals, 6 I 2 and 7(5, 0) = j(s). In 
what follows, we will denote j(s,e) = 7 e (s). 

Let W e be the vector field along 7 £ given by 



WJs) 



fryJs) 



dr 



s,e) 



Note that the vector fields W e and % commute 



d'y c?7 





' d d~ 




de' ds 



0. 



A variation r y £ of a horizontal curve 7 is called admissible if all curves 
7 e : I\ — > M are horizontal, 7 e (a) = 7(a) and 7 e (6) = 7(6) for all £ e J 2 - 
Observe that for an admissible variation of 7, the vector field Wq vanishes 
at the endpoints of 7: W (7(a)) = W (7(6)) = 0. 

Let us study an admissible variation 7 £ of a horizontal curve 7 in the case 
of H 1 , with the Riemannian metric defined in the previous Subsection. Since 
the variation is admissible, we have 

(%,Z x ) = (%,Zj) = (%,Z K ) = 0. 

In what follows, for an arbitrary vector field X on H 1 , we will denote by 
Xh and Xy the orthogonal projections of X to the horizontal distribution 
V C TH 1 and the vertical bundle span{Zj, Zj, Zic} respectively. 

The horizontality conditions (%, Z r ) = 0, for r = 1,J ', /C, yield 



de 



(%, Z r ) = (V Wo i, Z r ) + (7, V Wo Z r ) 



e=0 



= (V^W ,Z r ) + (j,V WoH Z r ) 

= 7 (Wo, Z r ) - (Wo, y-yZ r ) + (7, J r (WoH)) 

= j{W , Z r ) - (W 0h , J r m - ( J P (7), W 0fl ) 

= 7<W Q ,Z r )-2<Wo H ,J r (7)>, 

where we have used equation (|T5|) and V^ s Z r = 0. 
In fact the converse statement also holds. 

Lemma 1. Let W be any C 1 vector field along 7 such that ^(7(0)) 
W(7(6)) = and that satisfies 

= i(W,Z r )-2(W H ,Mj)). 
Then there exists an admissible variation 7 £ of 7 such that 

d 



de 



W. 



e=0 
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Proof. Let us decompose W = fj + W, with W-Lj and /(7(a)) = f(j(b)) = 
0. With this definition, we have 

(W,<y) = f, (W,J r m = (W,J r (j)), (W,Z r ) = (W,Z r ). 

Observe that the term fj will not contribute to any admissible variation, 
therefore we can assume that 1VJ_7. Let s e I\ and e > sufficiently small. 
Define the mapping 

F(s,s) = exp l{s) {eW{s)), 

where exp is the exponential map associated to the metric (-, •) of H 1 . 

If W is horizontal in some nonempty interval I C ii, then W = Wh and 
also (Wh, J r (j)) — ^i(Wn, Z r ) = 0. This implies Wh = A(p)7, but since 
W H ±i, then W H = 0. 

If W(sq) is not horizontal, then F(s,e) defines locally a surface which is 
foliated by horizontal curves and it is transversal to the horizontal distribu- 
tion, since it contains curves in nonhorizontal directions. This implies there 
exists a C 2 function g(s,e) such that 



7 B (s) = exp j(s) (g(s,e)W(s)) 

d 



f(so,e) = 1, we obtain 



e=0 



is a horizontal curve. Choosing q such that 

de 

an admissible variation 7 e of 7 with associated vector field W . □ 

With this result at hand, we can formulate the main theorem of this 
section. 

Theorem 3. Let 7 : [a, b] — > H 1 be a horizontal curve, parameterized by arc 
length. Then 7 is a critical point of the length functional if and only if there 
exist Aj, Xj, Xjc G R satisfying the second order differential equation 

(16) V^ 7 -2 A r J r (7)=0. 

r=X,J,K. 

Proof. Let 7 : I — [a, b] — > H 1 be a horizontal curve, parameterized by arc 
length, and let j e be an admissible variation of 7, with vector field U. The 
first variation of the length functional, see [6], is given by 



Hie) = ~ /< 
=0 Ji 



Suppose 7 is a critical point of the first variation, that is 

(Vyy, u) = 0. 
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The condition ||7|| = 1 implies (Vyy, 7) = 0. Since {j,Z r ) = for 
r = X,J ', ]C, we have 

= 7( 7 ,Z r ) = (V^Z r ) + (7,V^) 
= (V<vt, Z r ) + (7, J r ( 7 )} 
= (V^7, Z P ). 

Therefore, counting dimensions 

(18) V^ 7 = 9r(l)Jr(i). 

r=X,J,K 

In order to prove that the functions g r are constant, fix three C 1 functions 
fx, fjJz-.I^R such that / r (o) = f r (b) = and J f r = for r = X, J, /C. 

Consider a vector field C7 such that Ujj = Yl r =xjK frJr{j) an d (U, Z r )(s) = 
2 J s f r (t)dt. We claim that [/ satisfies 

j(U,Z r ) = 2{U H ,J r m, 
for r = X, ,7, /C. To see this, observe that 

7<M r ) = j- s UJ' f r (t)dt \ = 2f r (s) 

and also 

2(U H ,J r (j)) = 2/ £ f s J s (j),J r (j)\=2f r (s). 

\a=X,J,K I 

Thus, by Lemma [JJ we can conclude that U is a vector field for an admis- 
sible variation of 7. By the variational identity ( ITT)) , we obtain the equality 

= I (Vyy, U)= J2 f fr(Vji, Jr(j)), 

J* r=X,J,K^ 1 

which is valid for any three functions with mean zero. This implies that the 
functions (Vyy, J r (i)) are constant, and thus we obtain equation ([TBI) . for 
suitable constants Aj, Xj, X;c € K.- 

Conversely, let us assume that 7 is a horizontal curve, such that ||7|| = 1 
and it satisfies the differential equation (TT6|) . for some \x,\j,\/c G We 
need to show that 

for any C 1 -smooth vector field U, vanishing at the endpoints of 7 and sa- 
tisfying 

j{U,Z r ) = 2{U H ,J r (ff)), 

where r = X,J fC. 
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Let us write U 
9r{l{b)) = 0, then 



7. The intrinsic sub-Laplacian for S 7 with growth vector (6, 1) 

In pQ the authors presented an intrinsic form of the sub-Laplacian, by 
means of Popp's measure fi sR , introduced in [15]. The aim of this section is 
to construct this differential operator for the case of S 7 endowed with the 
contact distribution, introduced in Section 3. 

7.1. Construction of the intrinsic sub-Laplacian. Let (M,H,(-,-) s r) 
be a sub-Riemannian manifold, where H is a regular distribution. The basic 
idea is to define the intrinsic sub-Laplacian A sR f of a function / : M — >■ R 
of class C 2 , in analogy to the Riemannian case. To do this, let us define the 
horizontal gradient V sR f by the equation 

(19) (V sR f(p),v) sR = d p f(v), 

and the sub-Riemannian divergence div sR X of a horizontal vector field X 
by 

(20) div sR Xn sR = LxHsH, 

where fj, sR e f\ n (T*M) is a fixed non- vanishing n— form, known as Popp's 
volume form, and Lx denotes the Lie derivative in the direction of X. The 
intrinsic sub-Laplacian is given by 

(21) A sR f = dw sR (V sR f). 
For full details about its construction, see [H [15] . 

Remark: In the Riemannian case this definition coincides with the classical 
definition of the Laplacian, see for example [18]. As pointed out in [I], the 
regularity hypothesis over the distribution cannot be avoided since for exam- 
ple, in the case of the Grushin plane, the operator (1211) is not hypoelliptic. 

Let {Xi, . . . , Xk] be a local orthonormal basis of TL C TM and consider 
the corresponding dual basis {dXi, . . . ,dX k }. It is possible to find vector 



= U H + U v = U H + 22 9rZ r , where g r (l(a)) 

r=X,J,K 

-2 K f(Mi),u) = -2 Yl x r f(Mi)>u H ) 

r=X,J,K ^ I r=X,J,K " 1 

Kfj(u,z r ) = - Kfi(u v ,z r ) 

r=X,J,K ^ r=X,J,K ° 1 

r=X,J,K Jl r=X,J,K Ja 



□ 
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fields {X k+ i, . . . , X n } such that span{Xi, . . . , X n } = TM and such that 
Popp's volume form is locally given by 

(22) fi sR = dXi A ... A dX k A dX k+1 A ... A dX n . 

In this setting, the sub-Laplacian A sR f can be written explicitly as 

k / n 

(23) A sR f = J2 [L 2 X J + L Xr fJ2dX s ([X r ,X s \) 

r=l \ s=l 

7.2. Examples. The case of the intrinsic sub-Laplacian for S 3 is implied by 
the following result, characterizing Popp's volume form for contact manifolds 
of dimension 3. 

Proposition 9 ([TJ [15]). Let M be a three dimensional orientable contact 
manifold with a sub-Riemannian metric defined on its contact distribution. 
Let {Xi, X2} a local orthonormal frame for its contact distribution. Let 
X 3 = [Xi,X 2 ] and {dXi,dX 2 ,dX 3 } be the dual basis to {Xx, X 2 , X 3 } . Then 
the form dXi A dX 2 A dX 3 is an intrinsic volume form. 

In particular, for the sphere S 3 endowed with the contact distribution 
generated by the globally defined vector fields (jHJ), with commutator 

[X, Y](x) = 2V(x) = 2(-x 1 d X0 + x d xi - x 3 d X2 + x 2 d X3 ), 

Popp's volume form, as constructed above, is 2dX A dY A dV, and the 
intrinsic sub-Laplacian is given by 

A sR f = (X 2 + Y 2 )f. 

In general, we can extend the previous result to construct locally Popp's 
volume form over contact manifolds of arbitrary dimension. Let M be a 
contact manifold of dimension 2n+ 1, with contact form oj and contact dis- 
tribution £ = kerw. The distribution £ is bracket generating of step two, 
see [H]. Assume that M has a Riemannian metric g such that, in a neighbor- 
hood of each p e M, there is an orthonormal basis B = {vi, . . . , v 2n , v 2n+ \} 
for T P M satisfying £ p = span{f 1, . . . , v 2n }. Following the construction in [15] 
we have that Popp's volume form in this case is given locally by 

(24) fj, aR = 7Ti A . . . A7T 2n+ i, 

where B* = {iri, . . . , 7r 2n+ i} is the dual basis for B. 

In the case of the contact structure of S 7 , let us consider the vector fields 
Xx, . . . ,Xf presented in the Appendix. Since the vector fields X± and V4 
from equation (JS]) coincide, the contact distribution on S 7 introduced in 
Section 3 corresponds to 

H = keru = span{A 2 , . . . , X 7 }. 

In this context we have the following 
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Theorem 4. Let % be the contact distribution for S 7 and (•, -) s r the re- 
striction of the usual Riemannian metric in IR 8 to H. Then the intrinsic 
sub-Laplacian of (S 7 ,!!, (•, is given by the sum of squares 

a=2 

Proof. The construction of Popp's measure leads to the globally defined n 
form 

HsR = dX\ A ... A dX-?, 

which is precisely the Riemannian volume form of S 7 . Simple calculations 
show that 

(25) dX b ([X a ,X b ]) = (X b ,[X a ,X b ]) sR = 0, a = 2,..., 7 6=1,... 7. 
The theorem follows from formula (123]) . □ 

Remark: A complete list of the commutators [Jf ,Jff,], for a < b, can be 
found in (9j Section 8]. This list can be used to check equation ( 1251) directly. 

8. Heat operator for S 7 with growth vector (6, 1) 

The aim of this section is to show that the above constructed operator 
A s r commutes with the operator X\ . A similar observation was exploited 
to study the heat operator for the sub- Riemannian structure of SU(2) = S 3 

The main result of this Section is formulated as follows. 
Theorem 5. The operators A sR and X\ commute. 
Proof. Let us introduce the following change of coordinates for S 7 : 

xq + ix\ = e^ 1 cos r/i cos if) 



x 2 + ix 3 = e l & sin r/i cos ip 
X4 + ix$ = e^ 3 cos r]2 sin ip 
Xq + ix-j = e^ 4 sin r] 2 sin %p 

By the chain rule, the symbol of the sub-Laplacian A s ^ = X\ + . . . + Xj 
is a quadratic form with matrix 

//ii(77i,V) -1 -1 -1 \ 

-l h 2 (m,i } ) -i -looo 

-1 -1 h 3 {ri2,il>) -1 

-1 -1 -1 h 4 (r] 2 ,ij) 

sec 2 ^ 

csc 2 ^ 

\0 001/ 

where the coefficient functions hi,h 2 , and are given by 
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hi{vi,i>) = _ sec2 fa) Sec2 W ( _ 6 + 2 o ()K (2/ /j ) + a,s(2(/ /j - r)i - 



+2 cos(2^) + cos(2(?7i + ^)) 



-2 cos(2^) + cos(2(?7i + ^)) 

= Sec2 (^) CSc2 ^) ( 6 - 2 C0S(2 % ) + COSpfo - ^)) + 



+2 cos(2^) + cos(2(^ 2 + if>)) 

h 4 ( m , = _ CSC2 (^) CSC 2 W / _ g _ 2 CQs(2r?2) + cog(2(r?2 _ 



-2cos(2^) + cos(2(r7 2 + ^))J. 

Observe that hi,...,h± are independent of £i, . . . , £4. On the other hand, 
the vector field X 1; written in the new coordinates, becomes 

Since the coefficients of A sR are independent of the variables £15 £2; £3 
and £4, it is clear that the operators A sR and X\ commute. The Theorem 
follows. □ 

Let us denote by e~ tAaR the semigroup of operators acting on L 2 y with 
infinitesimal generator A sR . The operator e~ tAsR is known as the sub- 
Riemannian heat operator. As a consequence of Theorem we get the 
announced result. 

Corollary 3. Denoting by A57 the Laplace- Beltrami operator in S 7 with 
respect to the usual Riemannian structure, we have that 

e -tA g7 = e -t(A sR +X'{) _ e -tA 3 R e -tXf 

Proof. Since A57 = A sR +Xf, we have by the commutativity of the operators 

(27) e~ tA s 7 = e "* (Asfl+x i 2) = e~ tA " R e~ tx * , 

yielding to the stated result. □ 

The theory of unbounded operators allows us to rephrase the result in 
Corollary [3] as: 

Corollary 4. The sub- Riemannian heat operator e~ tAsR is given by 
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9. Appendix: Tangent vector fields to S' 



Octonion multiplication induces the following orthonormal basis of TS 
with respect to the restriction of the inner product (•, •) from M 8 to the 
tangent space T P S 7 at each p e S 7 . 



X x {x) 


= ~xid xo 


+ x d Xl 


- x 3 d X2 


+ x 2 d xa 


x$d x ^ 


+ x 4 d X5 


- x 7 d X6 


+ x 6 d X7 


X 2 (x) 


= -x 2 d X0 


+ x 3 d Xl 


+ x d X2 


- xid X3 


XqO X4 


+ x 7 d X5 


+ x 4 d X6 


- x 5 d X7 


X 3 (x) 


x 3 d X0 


- x 2 d xi 


+ x\d X2 


+ x d X3 


+ x 7 d X4 


-\- XqO X5 


x$d X g 


x 4 d X7 


X 4 {x) 


= x 4 d X0 


+ x 5 d xi 


-\- XqO X2 


- x r d X3 


+ x d X4 


- x x d X5 


~~ x 2 d Xf} 


+ x 3 d X7 


X 5 {x) 


= ~x 5 d X0 


x 4 d Xl 


- x 7 d X2 


xq9 X3 


+ x x d XA 


+ x d X5 


+ x 3 d X6 


+ x 2 d X7 


X 6 {x) 


x§d Xo 


+ x 7 d xi 


XaO X2 


+ x 5 d X3 


+ x 2 d X4 


- x 3 d X5 


+ x d X6 


- x x d X7 


X 7 (x) 


= -x 7 d X0 


- x 6 d Xl 


+ x 5 d X2 


+ x 4 d X3 


- x 3 d X4 


- x 2 d X5 


+ xid X6 


+ x d X7 
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